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We present a two-flavour linear sigma model with global chiral symmetry and vector and axial- 
vector mesons. We calculate irn scattering lengths and the decay widths of scalar, vector, and 
axial-vector mesons. It is demonstrated that vector and axial-vector meson degrees of freedom 
play an important role in these low-energy processes and that a reasonable theoretical description 
requires globally chirally invariant terms other than the vector meson mass term. An important 
question for meson vacuum phenomenology is the quark content of the physical scalar /o(600) and 
ao (980) mesons. We investigate this question by assigning the quark-antiquark a and ao states of 
our model with these physical mesons. We show via a detailed comparison with experimental data 
that this scenario can describe all vacuum properties studied here except for the decay width of the 
a, which turns out to be too small. We also study the alternative assignment /o(1370) and ao(1450) 
for the scalar mesons. In this case the decay width agrees with the experimental value, but the irir 
scattering length afj is too small. This indicates the necessity to extend our model by additional 
scalar degrees of freedom. 
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I. INTRODUCTION 



The fundamental theory of strong interactions, Quantum Chromodynamics (QCD), possesses an exact SU(3) C 
local gauge symmetry (the color symmetry) and an approximate global U(Nf)n x U(Nf)i, symmetry for Nf massless 
" quark flavours (the chiral symmetry). For sufficiently low temperature and density quarks and gluons are confined 
, into colorless hadrons (i.e., SU(3) C invariant configurations). Thus, it is the chiral symmetry which predominantly 
■^j" ' determines hadronic interactions in the low-energy region. 

, Effective field theories which contain hadrons as degrees of freedom rather than quarks and gluons have been 
developed along two lines which differ in the way in which chiral symmetry is realized: linear [l| and non-linear 
0. In the non- linear realization, the so-called non-linear sigma model, the scalar states are integrated out, leaving 
the pseudoscalar states as the only degrees of the freedom. On the other hand, in the linear representation of the 
symmetry, the so-called linear sigma model, both the scalar and pseudoscalar degrees of freedom are present. 

In this work, we consider the linear representation of chiral symmetry. An exactly linearly realized chiral symmetry 
implies that the QCD eigenstates come in degenerate pairs, the so-called chiral partners. Chiral partners have the 
same quantum numbers with the exception of parity and G-parity - for example, the scalar states sigma and pion 
and the vector states p and ai, respectively, are chiral partners. Experimental data in vacuum and at sufficiently 
low temperatures and densities of matter, however, show that the mass degeneracy is lifted, because the chiral 
U(Nf)n x U(Nf)L = U(i)v x U(1)a x SU(Nf)y x SU(N/)a symmetry is broken in two ways: explicitly and 
spontaneously. 

Due to the U(1)a anomaly [3|, the U(Nf)n x U(Nf)L symmetry is broken explicitly by quantum effects to U(l)v x 
SU(Nf)v x SU(Nf)A- In the case of small but nonzero degenerate quark masses, the latter is explicitly broken to 
U(Nf)y If the quark masses are not degenerate, the U{Nf)y symmetry is furthermore explicitly broken to U(l)v, 
corresponding to baryon number conservation. QCD also possesses discrete symmetries such as the charge conjugation 
(C), parity (P) and time reversal (X 1 ) symmetry (CPT), which are to a very good precision separately conserved by 
strong interactions. This fact offers further constraints in the construction of effective models of QCD. [A review of 
a possible, although small, CP violation in strong interactions may be found e.g. in Ref. Q.] 

In addition to the explicit breaking of axial symmetry SU(Nj)a due to nonzero quark masses, the latter symmetry 
is also spontaneously broken in vacuum by the non-vanishing expectation value of the quark condensate: (qq) = 
(QB,qL + q~L<lR) 7^ 5]. This symmetry breaking mechanism leads to the emergence of Nj — 1 pseudoscalar Goldstone 
bosons, as well as of massive scalar states representing the chiral partners of the Goldstone bosons. For Nf = 2, the 
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three lightest mesonic states, the pions, are identified with these Goldstone bosons of QCD. Their non-vanishing mass 
arises due to the explicit breaking of the chiral symmetry, rendering them pseudo-Goldstone bosons. 

In this paper we study an Nf = 2 linear sigma model which contains scalar (a, So) and pseudoscalar (77, 7?), and in 
addition also vector (cj, p) and axial-vector d\) degrees of freedom. Usually, such models are constructed under 
the requirement of local chiral invariance U(N /)# x U (JV/)^, with the exception of the vector meson mass term which 
renders the local symmetry a global one [a, |7[ ■ In a slight abuse of terminology, we will refer to these models as locally 
chirally invariant models in the following. A study of the QCD phase transition and its critical temperature T c within 
such a model can be found e.g. in Ref. |8j . 

However, as shown in Refs. [1, 0, the locally invariant linear sigma model fails to simultaneously describe 

meson decay widths and pion-pion scattering lengths in vacuum. As outlined in Ref. [T(| , there are at least two ways 
to solve this issue. One way is to utilize a model in which the (up to the vector meson mass term) local invariance 
of the theory is retained while higher-order terms are added to the Lagrangian [1, 0, ■ The second way which is 
pursued here is the following: we construct a linear sigma model with global chiral invariance containing all terms up 
to naive scaling dimension four [l2j |. The global invariance allows for additional terms to appear in our Lagrangian 
in comparison to the locally invariant case presented e.g. in Ref. Q. We remark that, introducing a dilaton field, one 
can argue fl3l . [l4| that chirally invariant terms of higher order than scaling dimension four should be absent. 

In Ref. we have presented a first study of meson decays and pion-pion scattering lengths in vacuum in the 
framework of the globally invariant linear sigma model. We have distinguished two different assignments for the 
scalar fields a = -^(uu + dd) and = -^(m — dd): (i) they may be identified with /o(600) and ao(980) which 
are members of a nonet that in addition consists of /o(980) and k(800); (ii) they may be identified with /o(1370) 
and ao(1450) which are members of a decuplet that in addition consists of /o(1500), ,fo(1710), and i^o(1430), where 
the additional scalar- isoscalar state emerges from the admixture of a glueball field [15]. In the following, we will 
refer to assignment (i) as Scenario I, and to assignment (ii) as Scenario II. In the latter, scalar mesons below 1 
GeV are not (predominantly) quark-antiquark states. Their spectroscopic wave functions might contain a dominant 
tetraquark or mesonic molecular contribution [16]. The correct assignment of the scalar quark-antiquark fields of the 
model to physical resonances is not only important as a contribution to the ongoing debate about the nature of these 
resonances, but it is also vital for a study of the properties of hadrons at nonzero temperature and density, where the 
chiral partner of the pion plays a crucial role jl7| . 



It is important to stress that the theoretical a and ao fields entering the linear sigma model describe pure quark- 
antiquark states, just as all the other fields (77, 7?, u>, p, fx, di). This property can be easily proven by using 
well-known large- N c results [l8j : the mass and the decay widths of both a and ao fields scale in the model as iV° and 
iV" 1 , respectively. 

In this paper we first investigate the consequences of Scenario I on various decay widths and pion-pion scattering 
lengths. This assignment is disfavored because a consistent description of all experimental data cannot be achieved. 
To reach this conclusion, vector and axial- vector degrees of freedoms play an important role. On the one hand their 
decays (such as p — > mr and a\ — > 777) and the role of the p meson in 7T7t scattering provide strong constraints, on 
the other hand they affect, indirectly but sizably, some decay channels, such as a — ¥ 7T7t. We then present a study of 
Scenario II. Although the latter is not yet conclusive because additional scalar fields (glueball, tetraquark) are not yet 
taken into account, our preliminary results for the decay widths (albeit not for the scattering length a®) are consistent 
with the data. 

The paper is organized as follows: in Sec. II we present the Lagrangian of our model and discuss the parameters 
which are known to very good precision and thus do not enter the fit of the decay widths and the scattering lengths. 
In Sec. Ill we present the formulas for the decay widths and the pion-pion scattering lengths which will be used to 
fit the remaining parameters and to compare the results to experimental data. This fit and comparison are discussed 
in Sec. IV, both for Scenario I and Scenario II. In Sec. V we summarize our results in the conclusions and give an 
outlook to future work. In the Appendix, we show the explicit form of our Lagrangian in terms of the meson fields. 
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II. THE LINEAR SIGMA MODEL WITH GLOBAL CHIRAL SYMMETRY 

A. The Lagrangian 

The Lagrangian of the globally invariant linear sigma model with U(2)r x U(2)l symmetry for Nf — 2 reads 

mil mug: 

C = Tr[(D Al $) t - moTr($ t $) - A 1 [Tr($ t $)] 2 - A 2 Tr($ t $) 2 

- jTt[(L^) 2 + (i?H 2 ] + ^Tr[(L^) 2 + (iF) 2 ] + Tr[iJ(* + $ f )] 

+ c(det $ + det $ f ) - 2 J5 2 (Tr {L^L", L"}} + Tr{R^[R", R"]}) 

- 2g 3 [Tr ({d^L v - teA^t 3 , L v \ + d v L^ - ieA„[t 3 , L„] } {U\ L v }) 

+ Tr {{d^R u - ieA^t 3 ,R u }+d^ - ieA y [t 3 , R^]} {R»,R U })] 

+ yTr($ t $)Tr[( J L Al ) 2 + (iF) 2 ] + /i 2 Tr[($iF) 2 + (L^$) 2 ] + 2h 3 Tr($R p & L>*). 

+ 54 {Tr [L^L U L„L U ] + Tr [R» R v R^R V ]} + g 5 {Tr [L^L tl L u L v ] + Tr [R» R„R U R u ] } 

+ g 6 Tr [R^Rfj] Tr [L V L V ] + g 7 {TrfL^L^] Ti[L v L v ] + Tv^R^] Tv[R"R„}} . (1) 

Note that the locally chirally invariant linear sigma model emerges from the globally invariant Lagrangian ([1} by 
setting hi=h 2 =h z = g 3 = 0, g 2 = g± = 9b = 96 = 9i = 9- 
In Eq. ©, 

$ = (cr + ii] N ) t° + (a + in) ■ t (2) 

contains scalar and pseudoscalar mesons, where t°, t are the generators of U (2) in the fundamental representation and 
T]n denotes the non-strange content of the 77 meson. Vector and axial- vector mesons are contained in the left-handed 
and right-handed vector fields: 

L» = (u* + fl l ) t° + + 3? ) ■ t , (3a) 
R» = K - /f ) t° + - a? ) • f , (3b) 

respectively. The covariant derivative 

= 9 m$ _ i gi (L^^ - $i?^) - ieA»[t 3 , $] (4) 

couples scalar and pseudoscalar degrees of freedom to vector and axial- vector ones as well as to the electromagnetic 
field A^. Note that local chiral invariance requires g\ = g. The left-handed and right-handed field strength tensors, 

L" v = d fl L u -ieA f *[t 3 ,L u }- {d u L" - ieA u [t 3 , L^}} , (5a) 
Ri iv = g » R » _ ieA ^ iR ^ _ {orgM _ ieA v [t 3 ,R' 1 ]} , (5b) 

respectively, couple vector and axial- vector mesons to the electromagnetic field A^ L . Explicit breaking of the global 
symmetry is described by the term Tr[7f($ + = ho<J (ho — const.). The chiral anomaly is described by the term 
c(det$ + det^) H. The model has been extended to include the nucleon field and its putative chiral partner; for 
details, see Refs. fl3l [20j. 

In the pseudoscalar and (axial-)vector sectors the identification of mesons with particles listed in Ref. [2l[ is 
straightforward, as already indicated in Eqs. ([2]) and (13111) - (| 3b I) : the fields 7? and ijn correspond to the pion and the 
SU(2) counterpart of the rj meson, r/N = (uu + dd)/V2, with a mass of about 700 MeV. This value can be obtained 
by "unmixing" the physical rj and rf mesons, which also contain ss contributions. The fields w M and p 5 " represent 
the w(782) and p(770) vector mesons, respectively, while the fields and di^ represent the /i(1285) and ai(1260) 
axial-vector mesons, respectively. (In principle, the physical tu and f± states also contain ss contributions, however 
their admixture is negligibly small.) Unfortunately, the identification of the a and do fields is controversial, the 
possibilities being the pairs {/o(600), a (980)} and {/o(1370), ao(1450)}. As mentioned in the Introduction, we will 
refer to these two assignments as Scenarios I and II, respectively. We discuss the implications of these two scenarios 
in the following. 

One may raise the question whether vector meson dominance (VMD) is still respected in the globally invariant 
linear sigma model (UJ. As outlined in Ref. [22j . there are two ways to realize VMD in a linear sigma model. The 
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standard version of VMD was introduced by Sakurai [23j and considers vector mesons as Yang-Mills gauge fields [24| . 
The gauge symmetry is explicitly broken by the vector meson masses. Another realisation of VMD was first explored 
by Lurie [25| whose theory contained a Lagrangian which was globally invariant. It is interesting to note that Lurie's 
Lagrangian contained direct couplings of the photon to pions and p mesons, as well as a p-ir coupling. It was shown 
in R,ef. [221 ] that the two representations of VMD are equivalent if the p-ir coupling g equals the photon-p coupling 
g p (the so-called "universal limit"). It was also shown that, if the underlying theory is globally invariant, the pion 
form factor at threshold F 7T (q 2 = 0) = 1 for any value of the above mentioned couplings. On the other hand, in 

Sakurai's theory F 1T (q 2 = 0) ^ 1 unless one demands g p7T7! — g p , or other parameters are adjusted in such a way 
that F- K (q 2 = 0) = 1. In other words, for any globally invariant model, and thus also for ours, one has the liberty of 
choosing different values for the photon-p and p-n couplings, without violating VMD. 



B. Tree-Level Masses 



The Lagrangian (TTJ contains 16 parameters. However, the parameters <?& with k = 3, 7 are not relevant for the 
results presented here so that the number of undetermined parameters decreases to eleven: 



mo, Ai, A 2 , mi, gx, 52, c, h , hi, h 2 , h 3 



(6) 



The squared tree-level masses of the mesons in our model contain a contribution arising from spontaneous symmetry 
breaking, proportional to (ft 2 . The value (ft is the vacuum expectation value of the a field and coincides with the 
minimum of the potential that follows from Eq. (TTJ . The a field is the only field with the quantum numbers of the 
vacuum, J PC = ++ , i.e., the condensation of which does not lead to the breaking of parity, charge conjugation, and 
Lorentz invariance. The potential for the a field reads explicitly 



V(a) 



-(ml - c)<7 2 + i 



Ai + 



hnO 



and its minimum is determined by 

= 



dV_ 
ck7 



c+ Ai + 



(7) 



(8) 



Spontaneous symmetry breaking corresponds to the case when the potential V(<ft) assumes its minimum for a nonva- 
nishing value a ~ (ft 0. In order to determine the fluctuation of the a field around the new vacuum, one shifts it 
by its vacuum expectation value (ft ^ 0, a — > a + (ft. The shift leads also to n^-fi and if-di mixing terms and thus to 
non-diagonal elements in the scattering matrix. These terms are removed from the Lagrangian by shifting the /1 and 
Si fields as follows: 



/f ft + Zwd» V 



where we defined the quantities 



TV , «i 



9i4> 

<m 2 



+ Zwd^ir , i]n — > Zr/N 



Zt: 



Z:= 1- 



-1/2 



(9) 



(10) 



Note that the field renormalisation of t/n and 7? guarantees the canonical normalization of the kinetic terms. This 
is necessary in order to interpret the Fourier components of the properly normalized one-meson states as creation or 
annihilation operators [6|]. Note also that the p and w masses as well as the fi and ai masses are degenerate in the 
globally as well as in the locally invariant model. Once the shift a — > a + <ft and the transformations ([§]) have been 
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performed, the mass terms of the mesons in the Lagrangian ([1} read: 

A 2 * 



77i~ = 7775 — c + 3 I Ai + 



m. 



Z 2 



z 2 



ml + c + \ \i + 



A 5 



c+ Ai + 



ml + 2cZ 2 



'/1 



m p = m 1 + y(^i + h>2 + h 3 ) , 



m\ + glcf) 2 + — (hi + h 2 - h 3 ) 



(11) 
(12) 

(13) 
(14) 
(15) 
(16) 



In Appendix A we show the Lagrangian in the form when all shifts have been explicitly performed. From Eqs. (|15p 
and (|16|) we obtain: 

m 2 ai =m 2 p + g 2 4> 2 -h 3 4> 2 . (17) 
The pion decay constant, / T is determined from the axial current, 



The large- N c dependence of the parameters is given by 

5 i, ff2 cx N^ 2 
A 2 , h 2 , h 3 , c oc N' 1 
Ai, hi oc N~ 
777q, m 2 cx N® , 
h oc N}/ 2 



c ' 

2 



(18) 



(19) 



We remind the reader that a vertex of n quark-antiquark mesons scales as iV c n . As a consequence, the parameters 
gi , <?2 scale as N c , because they are associated with a three-point vertex of quark-antiquark vector fields (of the 
kind p 3 ). Similarly, the parameters A 2 , /i 2 , h 3 scale as -ZV" 1 , because they are associated with quartic terms such as 
7r 4 and ir 2 p 2 . The parameter c is suppressed by a factor N c although it enters quadratic mass-like terms. This is due 
to the fact that the axial anomaly is suppressed in the large- N c limit. As is evident from Eq. (|12l) . the 77 n meson 
would also be a Goldstone boson for N c — ¥ 00. The parameters Ai, hi also describe quartic interactions, but are 
further suppressed by a factor 1 /N c because of the trace structure of the corresponding terms in the Lagrangian. The 
quantities TOq, m 2 are mass terms and therefore scale as 2V°. Then the pion decay constant f n scales as Nc . The 
quantity ho scales as N^ 2 in order that scales as N® as expected. Note that without any assumptions about the 
cr, do, and /1, ai fields, we immediately obtain that their masses scale as N® and their decay widths as N^T 1 , as we 
shall see in the following section. Therefore, they must also correspond to quark-antiquark degrees of freedom. 

There are, however, also approaches to the phenomenology of low-lying axial- vector mesons, such as the one in Ref. 
[Ill, where the Bethe-Salpeter equation is used to unitarize the scattering of vector and pseudoscalar mesons. Here, 
the Bethe-Salpeter kernel is given by the lowest-order effective Lagrangian. This leads to the dynamical generation 
of resonances, one of which has a pole mass of 1011 MeV and is consequently assigned to the ai(1260) meson. This 
unitarized approach is used in Ref. [27] to study the large-iV c behaviour of the dynamically generated resonances, 
with the conclusion that the ai(1260) resonance is not a genuine quark-antiquark state. 

However, it was shown in Ref. [Tij that, while unitarizing the chiral Lagrangian by means of a Bethe-Salpeter study 
allows one to find poles in the complex plane and identify them with physical resonances, it does not necessarily allow 
one to make a conclusion about the structure of those resonances in the large- N c limit. In order to be able to draw 
correct conclusions, a Bethe-Salpeter study requires at least one additional term of higher order not included in the 
Lagrangian of Refs. [26|,|27|. Alternatively, the Inverse Amplitude Method of Refs. [28| can be used. 
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A very similar approach to the one in Refs. [26J, |27J was also used in Ref. [29J where a very good fit to the r decay 
data from the ALEPH collaboration [30| was obtained by fine-tuning the subtraction point of a loop diagram. Note, 
however, that detuning the subtraction point by 5% will spoil the agreement with experimental data. Alternately, 
these data may be described by approaches with the ai(1260) meson as an explicit degree of freedom, such as the 
one in Ref. [12J, where ai(1260) is a quark- antiquark state and where the experimental ai(1260) spectral function 
is fitted very well. In Ref. [H], m ai (i260) — 1150 MeV and a full width r Ql (i260) — 410 MeV are obtained. Note 
that our results, as will be shown later, give very good results on the ai(1260) phenomenology, for example in the 
ai (1260) — > 7T7 and ai(1260) —> pir decay channels, see Sec. IV. A. 3. 

For the following discussion, it is interesting to note that the p meson mass, 

<j) 2 

m p = m 1 + — {hi + h 2 + h 3 ) , 

can be split into two contributions: the term m 2 which does not depend on the chiral condensate, and the term 
2 (^l + ^2 + /13) which depends quadratically on the condensate and vanishes in the chirally restored phase. We shall 

require that none of the two contributions be negative: in fact, a negative m\ would imply that the system is unstable 

,2 

when cj) —¥ 0; a negative ^-(hi + hi + ha) would imply that spontaneous chiral symmetry breaking decreases the p 
mass. This is clearly unnatural because the breaking of chiral symmetry generates a sizable effective mass for the light 
quarks, which is expected to positively contribute to the meson masses. This positive contribution is a feature of all 
known models (such as the Nambu-Jona-Lasinio model and constituent quark approaches). Indeed, in an important 
class of hadronic models (see Ref. [3l| and refs. therein) the only and obviously positive contribution to the p mass is 
proportional to 4> 2 (i.e., m x = 0). 



C. Equivalent set of parameters 

Instead of the eleven parameters in Eq. ([6]). it is technically simpler to use the following, equivalent set of eleven 
parameters in the expressions for the physical quantities: 

m^, rn„, m ao , m VH , m p , m ai , Z, <p, g 2 , hi, h 2 . (20) 

The quantities m^^mp, m ai are taken as the mean values for the masses of the tt, p, and ai meson, respectively, as 
given by the PDG [21]: m v = 139.57 MeV, m p = 775.49 MeV, and m ai = 1230 MeV. While m^ and m p are measured 
to very good precision, this is not the case for m ai ■ The mass value given above is referred to as an " educated guess" 
by the PDG [21|. Therefore, we shall also consider a smaller value, as suggested e.g. by the results of Ref. [12]. We 
shall see that, although the overall picture remains qualitatively unchanged, the description of the decay width of ai 
into pir can be substantially improved. 

As outlined in Ref. [ill ] , the mass of the rjN meson can be calculated using the mixing of strange and non-strange 
contributions in the physical fields r\ and 77' (958): 

V = T]N cos <p + r] S sin ip, T)' = -n N sin <p + ?] S cos ip, (21) 

where rjs denotes a pure ss state and tp ~ —36° [32|. In this way, we obtain the value m VN — 716 MeV. Given the 
well-known uncertainty of the value of <p, one could also consider other values, e.g., p = —41.4°, as published by the 
KLOE Collaboration [33|. In this case, m VN = 755 MeV. The variation of the 77 mass does not change the results 
significantly. 

The quantities <j> and Z are linked to the pion decay constant as <j>/Z = fir = 92.4 MeV. Therefore, the following 
six quantities remain as free parameters: 

m a , m aa , Z, g 2 , hi, h%. (22) 

The masses m CT and m ao depend on the scenario adopted for the scalar mesons. 

At the end of this subsection we report three useful formulas which link the parameters gi, and mi of the 
original set (j6j) to the second set of parameters (|20| [see also Eq. ([T0|l ] : 



9x=9i{Z) = ^-\l--^, (23) 



m ai 
ZfJ 


V_L 


mi 


( m l 







m\ = m\(Z, hi,h 2 ) = \[m 2 p + - Z 2 f 2 (g\ + h x + h 2 )] . (25) 
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III. DECAY WIDTHS AND tttt SCATTERING LENGTHS 



In this section, we quote the formulas for the decay widths and the tttt scattering lengths and specify their dependence 
on the parameters m a , m ao , Z, g 2 , h\, and h 2 . Using the scaling behavior (|19p we obtain that all strong decays and 
scattering lengths scale as N^ 1 , as expected. 

For future use we introduce the momentum function 



k(m a , m b , m c ) = Jm*- 2m 2 a (ml + m 2 ) + (ml - m 2 c ) 2 9(m a — m,b — m c ). (26) 

In the decay process a — > b + c, with masses m a , m^, m c , respectively, the quantity k(m a , mi,, m c ) represents the 
modulus of the three- momentum of the outgoing particles b and c in the rest frame of the decaying particle a. The 
theta function ensures that the decay width vanishes below threshold. 



A. The p — »• tttt decay width 



The decay width for p — > tttt reads 

1 p— >7T7T 



m . 



48*r< 



( 2m w 
V m p 



3/2 



gi z 2 + (1 - z 2 ) 



2\ 92 



(27) 



The experimental value is r^vV = (149.1 ± 0.8) MeV [21| . The small experimental error can be neglected and the 
central value is used as a further constraint allowing us to fix the parameter g 2 as function of Z: 

2 



92 = 92(Z) 



Z 2 -l 



9 4m;? 
gi Z 2 ± 



3tiT 



(cxp) 



p— >7T7T 



m 



y (m 2 — Ami) 3 / 2 



(28) 



Note that all input values in Eq. (|28|) are experimentally known [21|. The parameter g\ = gi(Z) is fixed via Eq. (|23l) . 

As apparent from Eq. ([2Sj) , two solutions for g 2 are obtained. The solution with the positive sign in front of the 
square root may be neglected because it leads to unphysically large values for the a\ — > pi: decay width, which is 
another quantity predicted by our study that also depends on g 2 [see Eq. (I59"|) ]. For example, the value Z = 1.6 
(see below) would lead to g 2 = 40 which in turn would give r ai _j. p7r = 14 GeV - clearly an unphysically large value. 
Therefore, we will take the solution for g 2 with the negative sign in front of the square root. In this case, reasonable 
values for both g 2 (see Table HJ and r ai _j. p7r (see Sec. IV. A. 3) are obtained. 



B. The /i — >■ aoTT decay width 



The decay width /i — > aQir reads 



r/i->a 7r( m a ! h 2 ) 



g\Z 2 A: 3 (m/ 1 ,m ao ,m 7r ) 



2tt 



:(h 2 



(29) 



There is a subtle point to comment on here. When the quark-antiquark ao state of our model is identified as the 
ao(980) meson of the PDG compilation (Scenario I), then this decay width can be used to fix the parameter h 2 as 

function of Z, h 2 = h 2 (Z), by using the corresponding experimental value T^*J^ ao7T — (8.748 ± 2.097) MeV (2lj . 



h 2 (Z) 



i „ m f, m„ 

3 ,2 ± h a 



(cxp) 



gi Z \ k 3 (m fl 



(30) 



Again, there are two solutions, just as in the case of the parameter g 2 . How strongly the somewhat uncertain 
experimental value of Tf 1 ^ a0 T T influences the possible values of h 2 , depends on the choice of the sign in front of the 
square root in Eq. (130[) . Varying Tf 1 ^ fao7r within its experimental range of uncertainty changes the value of h 2 by an 
average of 25% if the negative sign is chosen, but the same variation of r f 1 ^ aa iv changes h 2 by an average of only 6% 
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if the positive sign is considered. This is due to the fact that the solution with the positive square root sign yields 
larger values of h 2 ~ 80, while the solution with the negative sign leads to h 2 ~ 20. The absolute change of h 2 is 
the same in both cases. Our calculations have shown that using the negative sign in front of the square root yields 
a too small value of the 77-77' mixing angle ip = —9°. This follows by inserting h 2 into Eq. p"2"|) so that it is removed 
as a degree of freedom (i.e., replaced by Z) and calculating the mixing angle ip from Eq. pip using the experimental 
value of the ao — > r/n decay amplitude from Ref. 34] . For this reason, we only use the positive sign in front of the 
square root in Eq. (I31?)) . i.e., the constraint leading to higher values of h 2 . Then ip = —41.8° is obtained, in very good 
agreement with the central value quoted by the KLOE collaboration (33|, <p = —41.4° (see also Section IV.A.l). 

It may be interesting to note that only the (disregarded) lower value of h 2 leads to the expected behaviour of the 
parameter hi which [according to Eq. (fH?j) ] should be large- N c suppressed: the lower value of h 2 yields hi = 1.8 
whereas the higher value of h 2 yields hi = —68 (see Table [J). 

Note that if the quark- ant iquark ao meson of our model is identified as the ao(1450) meson of the PDG compilation 
(Scenario II) then the described procedure of replacing h 2 by Z using Eq. (|30)) is no longer applicable because the 
decay fi — > aoir is kinematically not allowed and its counterpart ao — > fin has not been measured. 



C. The ao — >■ r)Tr and ao — >■ r/'n decay amplitudes 



Our Nf = 2 Lagrangian contains the unphysical field rj^. However, by making use of Eq. (|21l) and invoking the 
OZI rule, it is possible to calculate the decay amplitude for the physical process ao — ► r/n as 



A aQ r}-n — COS (p A, 



From Eq. (fl} the formula for the decay amplitude containing the non-strange 77 field is 



A„ 



f TT(m ao Z, h 2 ) 



2 mi 
"1 



(m 2 a 



(31) 



(32) 



Note that Eq. (|32[) contains the unmixed mass m 7lN which enters when expressing the coupling constants in terms of 
the parameters (|20|) . as well as the physical mass m v = 547.8 MeV. The latter arises because the derivative couplings 
in the Lagrangian lead to the appearance of scalar invariants formed from the four-momenta of the particles emerging 
from the decay, which can be expressed in terms of the physical (invariant) masses. 

The decay width r oo _>,j follows from Eq. (|3"Tj) by including a phase space factor: 



qA m a : Z t ^2) 



fc(m ao , 777^,777^) 



(m aa ,Z, h 2 )Y 



(33) 



In the case of Scenario I, in which ao = ao(980), we shall compare the decay amplitude A aoV7r , Eq. pip, with 

the corresponding experimental value deduced from Crystal Barrel data: A^^J — (3330 ± 150) MeV (34J- This is 
preferable to the use of the decay width quoted by the PDG [2l[ for ao(980), which refers to the mean peak width, 
an unreliable quantity due to the closeness of the kaon-kaon threshold. 

In the case of Scenario II, in which ao = ao(1450), it is also possible to calculate the decay width ao(1450) — > 77'^, 
using the OZI rule. The amplitude A ao7) ' 7r (m ao , Z,h 2 ) is obtained following the same steps as in the previous case, 
Eq. ([32): 



A ao rj'Tr(lTla , Z, h 2 ) — — 



sin ip 



2 2 

777,„ — 777„ 



1 



2 m ai 



< 2 2 „.2 \ 

(m -m w -m ,) 



(34) 



where the difference compared to Eqs. p~T|) and (|32]l is the prefactor — sin tp and the physical 77' mass m v > — 958 MeV. 
The corresponding decay width reads: 



^a (1450)-s.?)'7r( TO ao ! Z, h 2 ) 



k(m ao ,m vl ,m w ) 2 
[A aov >TT (m ao , Z, h 2 )\ 



8ir 777,2 



(35) 



D. The ai — ^ TVf decay width 



We obtain the following formula for the a x — > 7r7 decay width: 

„2 



-^01—^^7 \Zj 



96tt 



(Z'~l)m ai 



(36) 
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Note that the a\ — > 7r7 decay width depends only on the renormalisation constant Z. Using rJi^Vy = (0.640 ± 0.246) 
MeV [2lj], one obtains Z — 1.67 ± 0.2. Unfortunately the experimental error for the quantity r ai _ i . 7r7 is large. Given 
that almost all quantities of interest depend very strongly on Z, a better experimental knowledge of this decay would 
be useful to constrain Z. In the study of Scenario I this decay width will be part of a \ 2 analysis, but still represents 
the main constraint for Z. 



We obtain the following formula: 



E. The a — Y -Kit decay width 



mir(m a ,Z, hi,h 2 ) 



32irm n 



2m T , 



9lZ*U 



-(h 1 + h 2 + h 3 ) 



(m 2 a -2 ml) 



(37) 



It is apparent from Eqs. (fl9[) that the sigma decay width decreases as the number of colours N c increases. Thus, the 
sigma field in our model is a qq state [281 ] . In Scenario I we have assigned the a field as /o(600), correspondingly we 
are working with the assumption that /o(600) [as well as ao(980)] is a qq state. In Scenario II, the same assumption 
is valid for the /o(1370) and ao(1450) states. 

Note that in Eq. (|37[) the first term in braces arises from the scalar aim vertex, while the second term comes from 
the coupling of the a to the a\, which becomes a derivatively coupled pion after the shift (HJ). Because of the different 
signs, these two terms interfere destructively As the decay width of a light a meson into two pions can be very well 
reproduced in the linear sigma model without vector mesons (corresponding to the case g\ — > 0), this interference 
prevents obtaining a reasonable value for this decay width in the present model with vector mesons, see Sec. IV. A. 2. 
This problem does not occur for a heavy a meson, see Sec. IV. B. 3 and Ref. [351 ] . 



F. The a\ — ► <j-k decay width 



The formula for the decay width reads 

Fa,i->-crrt( m <Ti Z, hi, h 2 ) 



k^(m ai , m a , m^) 2^2 
-9\Z 



6irm 6 ai 



-1 2 



h 3 ) 



(38) 



G. The ai — >• p-K decay width 



Let P be the four-momentum of the a\ meson, K\ the four-momentum of the p meson and K2 the four-momentum 
of the pion. Then the following formula for the a\ — > pir decay width is obtained: 



ai—tpir 



(Z) 



k(m a ^m p ,m v ) 



(V) 



m 2 p m 2 ai 



"-p "-at 

where h^ v is the vertex following from the relevant part of the Lagrangian ((T|) that reads 



(39) 



m a 1 



and 



K x -K 2 



m2 ai ~ ™> 2 P - m l 



P ■ K 2 = m ai En = m ai \Jh 2 + m% 
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Thus, we have 



( 2 2 

h%, = Z^fl A{g 2 - h 3 ) 2 + ^ [m 4 ai +mt + m 4 p + mlm 2 p + m 2 ai (m 2 - 2m 2 p ) + 3(r< - m 2 p - ml)m ai E^\ 

I m a 1 



o 9i92(9i - ha) , 2 2 2 , o v \ 



0152 (Si - M 



-4(m 0l - m p - m„)m ai E p E n ] 



("V ~ m*)m ai E p - 2m ai m p [E 7I + — 



E„ 



(V^) 2 = ^ {(.9 2 - fc) a < + |f-[« " <) 2 « +"4 " 2m 2 .) + {ml + <) 



2 \ 4 



-4(m Ql - m p - m^)m ai E p E^} 



5152(3? - M 



{h lu ,P'>KZy = ( g i-h 3 yz*fimi 1 E' p . 



H. The tree-level scattering lengths 

The partial wave decomposition [36j j leads to the following formula for the s-wave 1 = pion-pion scattering length 
Oq (in units of m~ x \. 



2 ~4 7T 



ag(Z,m CT ,/n) = -^l2^^ 

47T \ 771' 



16 



[120? - 2(fti + h 2 ) - 14/i 3 ] 



2 ryl 1 "V 



2m^ +m 2 p - —(h 1 + h 2 + ha) 



5 i? 2 777 2 — 777 2 

8 T 2 

Z 2 777 2 — 777? 



20 



1 



4m?. — to?. 



2 7 2, "'tt 



9\Z 



m p — —{h\ + ba + ha) 



Z 2 m 2 r — T77 2 . 

24 



(40) 



We use the value 



(cxp) 



= 0.218 ± 0.020 in accordance with the 2003 and 2004 data from the NA48/2 collaboration 



An analogous calculation leads to the s-wave 1 = 2 pion-pion scattering length a 2 ,: 



Oq(Z, 777 CT ,/7l 



1 / Z mz — 777 



4?r 



7T 1 „2 7 4 " t, 7T 



4/ 2 



5^ 



m 2 - —(hi + h 2 + h 3 ) 



2 7 2j. "'tt 



.9^ 



m p - — (/ii + ft, 2 + M 



Z 2 777?. 777 2 



20 



(41) 



The experimental result for a 2 , from the NA48/2 collaboration is au (exp) = -0.0457 ± 0.0125 [37|. Note that the tttt 
scattering lengths were also studied away from threshold in Ref. |38| . in a model quite similar to ours. 

IV. STUDY OF DIFFERENT SCENARIOS FOR THE STRUCTURE OF SCALAR MESONS 



In this section we discuss two different interpretations of the scalar mesons. The following subsection describes 
the results obtained when /o(600) and ao(980) are interpreted as scalar quarkonia (Scenario I). Then we discuss the 
results obtained when / (1370) and a (1450) are interpreted as scalar quarkonia (Scenario II). 
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A. Scenario I: Light Scalar Quarkonia 

1. Fit procedure 

As a first step we utilize the central value of the experimental result Tp*$n = 149.1 MeV l2ll in order to express 
the parameter g 2 as a function of Z via Eq. (f2"5)l . Moreover, we fix the mass m O0 = 0.98 GeV [21j and we also use the 
central value r/ 1 _ i . Qo7r (Z, h 2 ) = 8.748 MeV to express h 2 as a function of Z. The results are practically unaffected by 
the 6% uncertainty in h 2 originating from the uncertainty in Y j^ aolT , see Eq. (|30[) . 

As a result, the set of free parameters in Eq. (|2"2l is further reduced to three parameters: 

Z, TO,, h t . (42) 

Note that in this scenario the field a is identified with the resonance /o(600), but the experimental uncertainty on its 
mass is so large that it does not allow us to fix m a . We therefore keep m a as a free parameter. 

We now determine the parameters Z , hi, and m CT using known data on the a\ — > 7r7 decay width (|36[) and on the rnr 
scattering lengths a® and Oq reported in Eqs. (14T))) and (|4"Tj) . This is a system of three equations with three variables 
and can be solved uniquely. We make use of the \ 2 method in order to determine not only the central values for our 
parameters but also their error intervals: 



V Ai decay / iG{0,2} \ Afl / 

The errors for the model parameters are calculated as the square roots of the diagonal elements of the inverted Hessian 
matrix obtained from ^(Z, m a , hi). The minimal value is obtained for \ 2 = 0, as expected given that the parameters 
are determined from a uniquely solvable system of equations. The values of the parameters are as follows: 

Z = 1.67 ± 0.2 , m a = (332 ± 456) MeV , hi = -68 ± 338 . (44) 

Clearly, the error intervals for m a and hi are very large. Fortunately, it is possible to constrain the hi error 
interval as follows. As described at the end of Sec. II. B, the p mass squared contains two contributions - the bare 
mass term m\ and the quark condensate contribution (~ (f> 2 ). The contribution of the quark condensate is special 
for the globally invariant sigma model; in the locally invariant model m p is always equal to mi [8J. Each of these 
contributions should have at most the value of 775.49 MeV (= m p ) because otherwise either the bare mass or the 
quark condensate contribution to the rho mass would be negative, which appears to be unphysical. A plot of the 
function mi — mi (Z, hi,h 2 {Z)), sec Eq. ([25]), for the central values of Z = 1.67 and T^^^ = 8.748 MeV is shown 
in Fig. [[] 

Note that varying the value of T f^^ agn within its experimental boundaries would only very slightly change hi 

by ±4 and this parameter is thus unaffected by the experimental error for T^^. If the value of mi was known 
exactly, then Eq. (123)) would allow us to constrain hi via Z. However, given that at this point we can only state that 
< mi < m p , for each Z one may consider all values of hi between two boundaries, one obtained from the condition 
mi(Z, hi, h 2 (Z)) = and another obtained from the condition mi(Z, hi, h 2 (Z)) = m p . For example, using the central 
value of Z = 1.67, we obtain —83 < hi < —32. The lower boundary follows from mi = m p and the upper boundary 
from mi = 0, see Fig. [T] Note that the central value hi = —68 from Eq. ([44]) corresponds to mi = 652 MeV. If the 
minimal value of Z = 1.47 is used, then hi = —112 is obtained from mi = m p and hi = —46 from mi = 0. Thus, 
— 112 < hi < —46 for Z = 1.47. Analogously, —64 < hi < —24 is obtained for the maximal value Z = 1.87. 

Clearly, each lower boundary for hi is equivalent to mi = m p and each upper boundary for hi is equivalent to 
mi = 0. Thus, in the following we will only state the values of Z and mi; hi can always be calculated using Eq. (I25[) . 
In this way, the dependence of our results on mi and thus on the origin of the p mass will be exhibited. 

The value of m CT can be constrained in a way similar to hi using the scattering length a{j; the scattering length a§ 
possesses a rather large error interval making it unsuitable to constrain m a . Figure [2] shows the different values for 
ag and dp depending on the choice of Z and mi. 

It is obvious that the value of a[j is only consistent with the NA48/2 value [37}, if m CT is in the interval [288, 477] 
MeV, i.e., m CT = 332^44 5 MeV. This value for m a follows if the parameters Z and mi are varied within the allowed 
boundaries. If we only consider the Oq curve that is obtained for the central values of Z and mi, a much more 
constrained value of m a = 3321^ MeV follows from Fig. [2] We will be working with the broader interval of m a . 
Even then, constraining mi to the interval [0, m p ], the error bars for m a are reduced by at least a factor of three in 
comparison to the result (14"4"1) following from the % 2 calculation. 
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i«j(MeV) 




FIG. 1: mi as function of hi, constrained at the central value of Z = 1.67. The black dot marks the position of 

central values hi = —68 and mi = 652 MeV. 




leV) 



FIG. 2: Scattering lengths ajj and a 2 , as function of m a [the shaded band corresponds to the NA48/2 value of a°; no 

error interval is shown for a 2 , due to the large interval size [37j]. 



We summarize our results for the parameters Z and m a : 

Z = 1.67 ± 0.2 , m a = 332t\f MeV . 

The central values of all parameters of the original set (j6|) are given in Table HI They follow from the x 2 fit (m a , hi), 
via decay width constraints (/12, gi), and from Eqs. (ITT1) - (|16p and (j2"5)) - (IM|) . The central values of Z, m CT , and /ii, 
Eq. (|44|) , have been used to calculate all other parameters. We neglect the errors, apart from those of mi, which in 
this scenario vary in a large range. 



Parameter 


m a 


hi 


h 2 


h 3 


gi 


92 


m 


mi 


Ai 


A 2 


c 


ho 


Value 


332 MeV 


-68 


80 


2.4 


6.4 


3.1 


210 MeV 


652+^ MeV 


-14 


33 


88744 MeV 2 


1 ■ 10 b MeV a 



TABLE I: Central values of parameters for Scenario I. 



Note that the values of a§ depend strongly on the choice of the parameters Z and m\. Whereas for the central 
values of Z and mi this scattering length is constant and has the value Oq = —0.0454, its value increases if Z and mi 
are considered at their respective boundaries, see Fig. [2j 
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The value of Z alone allows us to calculate certain decay widths in the model. For example, as a consistency check 
we obtain r ai _ i . Tr7 = 0.640+o'23i MeV which is in good agreement with the experimental result. Also, given that the 
o-o —> riNir decay amplitude only depends on Z, it is possible to calculate the value of this amplitude, Eq. (|32j) . For 
Z = 1.67, we obtain A ao ^ ?)w = 3939 MeV for the decay amplitude ao — > i]ir involving the physical rj field if the 77-77' 

mixing angle of <p = —36° 32] is taken. The Crystal Barrel data [13 read A^S^n = 3330 MeV and hence there 
is an approximate discrepancy of 20%. If the KLOE Collaboration [33| value of ip = —41.4° is considered, then the 
value of Aa^r,?; — 3373 MeV follows - in perfect agreement with the Crystal Barrel value. From this we conclude 
that this scenario prefers a relatively large value of the 77-77' mixing angle. In fact, if we use the Crystal Barrel value 
A^a^r/TT = 3330 MeV as input, we would predict ip = —41.8° for the central value of Z as well as ip — —42.3° and 
ip = —41.6° for the highest and lowest values of Z, respectively, i.e., ip = — 41.8°+q'5o. This is in excellent agreement 
with the KLOE collaboration result tp = -41.4° ± 0.5°. 



2. The decay a — > nn 



The sigma decay width T^^^ depends on all three parameters Z, mi (originally hi), and m CT . In Fig. [3] we show 
the dependence of this decay width on the sigma mass for fixed values of Z and mi, varying the latter within their 
respective boundaries. 
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FIG. 3: F0._j.7r7r as function of m CT for different values of Z and mi. The PDG f- il q uotes rfc 1 ^ = (600 — 1000) 
MeV; the results from the chiral perturbation theory suggest T a ^Kir = 544 MeV [39| and r o ._ ) . 7r7r = 510 MeV [40| . 



Generally, the values that we obtain are too small when compared to the PDG data [2l| and to other calculations 
of the sigma meson decay width, such as the one performed by Leutwyler et al. [39| who found r o -_ > 7r T r/2 = 272+5* 2 5 
MeV and Pelaez et al. [40| who found = (255 ± 16) MeV. The largest values for the decay width that we 

were able to obtain within our model are for the case when Z is as small as possible, Z = 1.47, and mi = 0, i.e., when 
the p mass is solely generated by the quark condensate. As seen above, for this case the scattering lengths allow a 
maximum value m a = 477 MeV, for which Ta^^^ = 145 MeV. In all other cases, the decay width is even smaller. 
However, as will be discussed in Sec. IV. A. 3, the case mi = leads to the unphysically small value r ai ^ a7r ~ and 
should therefore not be taken too seriously. As apparent from Fig. (5J excluding small values of mi would require 
smaller values for m CT in order to be consistent with the scattering lengths. According to Fig. [3J however, this in turn 
leads to even smaller values for the decay width. 

Hence, we conclude that the isoscalar meson in our model cannot be /o(600), thus excluding that this resonance is 
predominantly a qq state and the chiral partner of the pion. Then the interpretation of the isospin-one state ao(980) 
as a (predominantly) quarkonium state is also excluded. The only choice is to consider Scenario II, see Sec. IV. B, i.e., 
to interpret the scalar states above 1 GeV, /o(1370) and ao(1450), as being predominantly quarkonia. If the decay 
width of /o(1370) can be described by the model, this would be a very strong indication that these higher-lying states 
can be indeed interpreted as (predominantly) qq states. Note that very similar results about the nature of the light 
scalar mesons were also found using different approaches: from an analysis of the meson behaviour in the large-V c 
limit in Refs. [28| and [4l| as well as from lattice studies, such as those in Refs. 
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We remark that the cause for preventing a reasonable fit of the light sigma decay width is the interference term 
arising from the vector mesons in Eq. p7[) . In the unphysical case without vector meson degrees of freedom, a 
simultaneous fit of the decay width and the scattering lengths is possible [35j | . 

3. Decays of the ai meson 

We first consider the decay width r^-^. For a given m Ql , this decay width depends only on Z. The PDG quotes 
a rather large band of values, Ta^pn — (250 — 600) MeV. For m ai = 1230 MeV, our fit of meson properties yields 
Z = 1.67 ± 0.2. The ensuing region is shown as shaded area in Fig. 01 For m ai = 1230 MeV, T ai ^ p7T decreases from 
2.4 GeV to 353 MeV, if Z varies from 1.47 to 1.87. 

We also observe from Fig. [4] that the range of values for Z, which give values for T^^p^ consistent with the 
experimental error band, becomes larger if one considers smaller masses for the a\ meson. We have taken m ai = 1180 
MeV and m ai — 1130 MeV, the latter being similar to the values used in Refs. 12) and [43j. Repeating our calculations, 
we obtain a new range of possible values for Z, Z ~ 1.69±0.2 for m ai — 1180 MeV and Z ~ 1.71. ±0.2 for m ai = 1130 
MeV. For the respective central values of Z we then compute r™ Q ±p" 8 ° M ° V = 483 MeV (^ ai =ii80MeV = 1 69 ) 
and r™% p ^ 13 ° McV = 226 MeV (2" l »i =1130McV = 1.71), in good agreement with experimental data. All other results 
remain valid when m ai is decreased by about 100 MeV. Most notably, the /o(600) decay width remains too small. 




z 



1.6 1.7 1.8 



FIG. 4: T^^pTr for different values of m ai . The shaded area corresponds to the possible values of r ai ^ p7r as stated 

by the PDG. 



We also consider the a\ — > air decay width. Experimental data on this decay channel [2l| are inconclusive. The 
value T ai ^ av = 56 MeV is obtained for the central values of Z, mi, m a , and Y f^ aoV (which was used to constrain hi 
via Z). Taking the limit m\ — pulls the value of T ai ^, av down to practically zero, regardless whether Z = Z m \ n or 
Z = Zmax- This is an indication that the mi = limit, where m p is completely generated from the quark condensate, 
cannot be physical. Note that the case Z = Z max = 1.87 and mi = m p , i.e., where the quark condensate contribution 
to the p mass vanishes, leads to a rather large value of r ai ^ a7T , e.g., for the central value of m a = 332 MeV the 
value of r ai _j. CTW = 120 MeV follows. Interestingly, this picture persists even if lower values of m Ql are considered. 
Improving experimental data for this decay channel would allow us to further constrain our parameters. 

4- The case of isospin-exact scattering lengths. 

So far, the values of the scattering lengths used in our fit, a[] = 0.218±0.020 and ag = — 0.0457±0.0125 (37|, account 
for the small explicit breaking of isospin symmetry due to the difference of the up and down quark masses. However, 
in our model the isospin symmetry is exact. Thus, one should rather use the isospin-exact values a®^ = 0.244 ±0.020 
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2 (j\ i 1 

and a l ' — —0.0385 ± 0.0125 [44|. In this section we will briefly show that the conclusions reached so far remain 



qualitatively unchanged if the isospin-exact values for the scattering lengths are considered. 

Performing the \ 2 fit, Eq. (031), with r ai _ i . 7r7 , aj] ™ and o,q ^ as experimental input yields Z = 1.67±0.2 - unchanged 
in comparison with the previous case (Z is largely determined by r ai _j.^ 7 which is the same in both x 2 calculations), 
hi = — 116 ± 70, and = (284 ± 16) MeV. Note that in this case the errors are much smaller than previously. The 
reason is that the mean value of to ct is almost on top of the two-pion decay threshold and thus leads to an artificially 
small error band. For such small values of m CT the decay width T a ^ n7T is at least an order of magnitude smaller than 
the physical value, but even for values of m CT up to 500 MeV (not supported by our error analysis) the decay width 
never exceeds 150 MeV, see Fig. [3] 



B. Scenario II: Scalar Quarkonia Above 1 GeV 



1. General discussion 



A possible way to resolve the problem of the unphysically small two-pion decay width of the sigma meson is to 
identify the fields a and ao of the model with the resonances /o(1370) and ao(1450), respectively. Thus, the scalar 
quarkonium states are assigned to the energy region above 1 GeV. In the following we investigate the consequences 
of this assignment. However, the analysis cannot be conclusive for various reasons: 

(i) The glueball field is missing. Many studies find that its role in the mass region at about 1.5 GeV is crucial, 
since it mixes with the other scalar resonances. 

(ii) The light scalar mesons below 1 GeV, such as /o(600) and ao(980), are not included as elementary fields in our 
model. The question is if they can be dynamically generated from the pseudoscalar fields already present in our 
model by solving a Bethe-Salpeter equation. If not, they should be introduced as additional elementary fields 
from the very beginning [see also the discussion in Ref. (14|]. 

(iii) Due to absence of the resonance /o(600), the hit scattering length q,q cannot be correctly described at tree- level: 
whereas stays always within the experimental error band, a,g clearly requires a light scalar meson for a proper 
description of experimental data because a large value of m a drives this quantity to the Weinberg limit (~ 0.159) 
which is outside the experimental error band. 

Despite these drawbacks, we turn to a quantitative analysis of this scenario. 



2. Decay of the ao(1450) meson 

As in Scenario I, the parameter g 2 can be expressed as a function of Z by using the p — > irn decay width (|27[) . 
However, the parameter h 2 can no longer be fixed by the f± — > aQir decay width: the ao meson is now identified with 
the ao(1450) resonance listed in Ref. [21fl . with a central mass of m ao = 1474 MeV, and thus fi is too light to decay 
into ao and ir. One would be able to determine h 2 from the (energetically allowed) decay ao(1450) — > /i7r, but the 
corresponding decay width is not experimentally known. 

Instead of performing a global fit, it is more convenient to proceed step by step and calculate the parameters 
Z, hi, h 2 explicitly. We vary m a = Ti/ (i37o) within the experimentally known error band |2l| and check if our result 
for rj ( 1 3 7 o)_i. 7r7r is in agreement with experimental data. 

We first determine Z from ai — > irj, Eq. (f3"6")) , and obtain Z — 1.67 ±0.21. We then immediately conclude that the 
ai — > pit decay width, Eq. (|39[) . will remain the same as in Scenario I because this decay width depends on Z (which 
is virtually the same in both scenarios) and g 2 [which is fixed via r p _,. W7r , Eq. (|28l) . in both scenarios]. 

The parameter hi, being large- N c suppressed, will be set to zero in the present study. We then only have to 
determine the parameter h 2 . This is done by fitting the total decay width of the ao(1450) meson to its experimental 
value 

r ao (1450)(^, h 2 ) = T aa ^ vri + T aa ^ m y + T aa ^KK + ^ a a ^uir-K = ^ao^iHO) = (^65 ± 13) MeV. (45) 

Although kaons have not been included into the calculations, we can easily evaluate the decay into KK by using 
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flavour symmetry 



r 17 h\ 9 H m a , m K,m K ) 2 
1 a (1450)-s-A"MA — ^ 87rr7?2 [ A a KK \4, "2jJ i 



A ao KK(Z,h 2 ) = 



2ZU 



1-5— r-fa-fc 

ai 



(m 



(46) 
(47) 



The remaining, experimentally poorly known decay width ^ ao (lA50)^uinn can be calculated from the sequential decay 
ao — > W/O — > aj7T7r. Note that the first decay step requires the p to be slightly below its mass-shell, since m ao < rrip + m^. 
We denote the off-shell mass of the p meson by x. From the Lagrangian ([lj we obtain the following formula for the 
ao — > cup decay width: 



r (™\ k{ma ,m u ,x) f -,2 7 2 ( 2 



87rmL 



3 - 



{ml -x 2 -ml) 



2 \2n 



The full decay width r ao ( 14 5o)_>. W7r7r is then obtained from the following equation: 

/>oo 

r ao (1450)^w7r7r = / dxT ag ^ LJp (x) d p (x) , 

Jo 

where d p (x) is the mass distribution of the p meson, which is taken to be of relativistic Breit-Wigner form: 



d P (x) = N 



2-p(oxp) 



X Z T 



p— >7T7T 



(x 2 - m*)* + (xrp&y 



(x - 2m 7r ) , 



where T 



(exp) 
p— >ir-K 



149.1 MeV and m p = 775.49 MeV [2T(. (In general, one should use the theoretical quantity T 

^(cxp) 



(48) 



(49) 



p — y7T7T ) 



which is itself a function of x, instead of Tp-nf^, see for instance Ref. [45j and refs. therein. This is, however, 
numerically irrelevant in the following.) The normalization constant N is chosen such that 



dxdp(x) = 1 , 



(50) 



in agreement with the interpretation of dxd p (x) as the probability that the off-shell p meson has a mass between x 
and x + dx. 

Inserting Eqs. (f3"3"|) . (|3"5|) . (14"6"1) . and (14"51) into Eq. (|4"5)) , we can express h 2 as a function of Z, analogously 
to Eq. (|28[) where 172 was expressed as a function of Z. Similar to that case, we obtain two bands for h 2 , 
— 115 < /12 < —20 and —25 < h 2 < 10, the width of the bands corresponding to the uncertainty in determining 
Z, Z = 1.67 ± 0.21. Both bands for h 2 remain practically unchanged if the 5% experimental uncertainty of 

r^^gps is taken into account and thus we only use the mean value 265 MeV in the following. Since hi is 
assumed to be zero, Eq. (|25[) allows to express m\ as a function of Z, m\ = mi(Z,hi = 0,h 2 (Z)) (we neglect 
the experimental uncertainties of m p , m ai , and f^). The result is shown in Fig. [5J The first band of (lower) h 2 
values should be discarded because it leads to mi > m p . The second set of (higher) values leads to mi < m p only 
if the lower boundary for Z is 1.60 rather than 1.46. Thus, we shall use the set of larger h 2 values and take the 
constraint mi < m p into account by restricting the values for Z to the range Z = LGT^q' 2 ^- As can ^ e seen from Fig. 
this sets a lower boundary for the value of mi, mi > 580 MeV. Thus, in this scenario we obtain mi = 720]ti4o MeV. 

The values for the other parameters can be found in Table [TX| (only central values are shown with the exception of 
mi where the corresponding uncertainties are stated as well). 



Parameter 


hi 


h 2 


h 3 


9i 


92 


mo 


mi 


Ai 


A 2 


c 


h 


Value 





4.7 


2.4 


6.4 


3.1 


-811987 MeV 


720+^ MeV 


-3.6 


84 


88747 MeV^ 


1 ■ 10 b MeV^ 



TABLE II: Central values of the parameters for Scenario II. 
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Note that Ai <C A2, in agreement with the expectations from the large- N c limit, Eq. (fl"9|) . The value of mi = 720 
MeV is sizable and constitutes a dominant contribution to the p mass. This implies that non-quark contributions, for 
instance a gluon condensate, play a decisive role in the p mass generation. 

As a final step, we study the ratios r oo ( 1450 )^ r) ' 7r /r oo (x45o)^^7r and r oo(1450) ^. K y/r ao ( 1450 )^ w . Their experimental 
values read [2l| 

p(cxp) p(oxp) _ 

a (1450Wx = Q.35 ± 0.1 6; ao(1450)-fJCX ^ Q.88 ± 0.23. (51) 

-p(oxp) -p(cxp) v ' 

1 a (1450)->rj7r a o (1450)-»J?7r 

Using the central value Z = 1.67 and ip = —36° for the rj — rf mixing angle, we obtain r ao ( 1 4 5 o)_j.^' W /r ao ( 1 45o)_ >r;7r = 
1.0 and r oo( - :L450 j_ >if -jf/r a[1 (- 1 45Q)_ > . r;7r = 0.96. The latter is in very good agreement with the experiment, the former a fac- 
tor of two larger. Note, however, that according to Eqs. (|3~Tj) and the value of the ratio r ao ( 1 4 5 o)_j. r) / 7r /r ao ( 1 4 5 o)_j. r)7r 
is proportional to sin 2 ip/ cos 2 p. If a lower value of the angle is considered, e.g., ip — —30°, then we obtain 
rao(i450)^j)'7r/r ao (i45o)^r(7r = 0.58 for the central value of Z and the central value of r ao (!4 5 o) in Eq. (l4"5l) . Tak- 



ing Z = Z max and the upper boundary r^ 450) = 278 MeV results in r Qo(1450 )^^/ 7r /r ao ( 1 4 50 )^, ;7r = 0.48, i.e., 
agreement with the experimental value. Therefore, our results in this scenario favour a smaller value of <p than the 
one suggested by the KLOE Collaboration [33[. 

It is possible to calculate the decay width r ao ( 14 5 )_j. W7r7r using Eq. (jig)) . We have obtained a very small 
value r aD (i 4 5o)_j, W7r7r = 0.1 MeV. From Eq. ([33]) we obtain r ran n ^bdV-^tt = 89.5 MeV, such that the ratio 
r ao(i450)^7r7r/r a o(i450)^r ) 7r = 0.0012, in contrast to the results of Ref. |4ej. 



m 



3. Decay of the /o(1370) meson 

It is now possible to calculate the width for the /o(1370) —> hit decay using Eq. ([3T)) . The decay width depends on 
the /o(1370) mass, Z, hi, and hi which is expressed via Z using Eq. (|4"5"|) . The values of the latter three are listed in 
Table HU In Fig. [S] we show the decay width as a function of the mass of /o(1370). 

Assuming that the two-pion decay dominates the total decay width, we observe a good agreement with the experi- 
mental values if m^ ( 137 o) < 1380 MeV. Other contributions to the decay width are likely to reduce this upper bound 
on toj ( 1370 ) somewhat. Nevertheless, the correspondence with the experiment is a lot better in this scenario where 
we have identified /o(1370) rather than / (600) as the (predominantly) isoscalar qq state. Note that this result has 
been obtained using the decay width of the ao(1450) meson (in order to express h,2 via Z), which is also assumed to 
be a scalar qq state in this scenario. 

It is remarkable that vector mesons are crucial to obtain realistic values for the decay width of /o(1370): without 
vector mesons, the decay width is ~ 10 GeV and thus much too large. This is why Scenario II has not been considered 
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in the standard linear sigma model. 

650 




FIG. 6: Dependence of the /o(1370) decay width on mj (i37o). The experimental value of the mass is expected to be 
in the range 1200-1500 MeV, the width between 200 and 500 MeV [H|. 

The four-body decay /o(1370) — > 47r can also be studied. Similarly to the ao(1450) — > cop decay, we view /o(1370) — > 
Air as a sequential decay of the form /o(1370) — s- pp — > Att. The Lagrangian (TTJ leads to 



r/ (i3ro)^ P(0 (aJi,^) = ^ ^2 (/H + /12 + /13) z f v 



m 



fo 



1 + A ( m / - x l - x lf 



m- 



+ 



4m4 



(52) 



where x% and £2 are the off-shell masses of the p mesons. The decay width T is then given by 



/o(1370)^4tt 



dxi dx 2 T fo{1370) ^ pp (x 1 ,x 2 ) d p (xi) d p (x 2 ) , 



with r 



/o(1370)^pp 



(21,2:2) from Eq. and d p (x) from Eq. ([35]). 



Using the previous values for the parameters we obtain that the pp contribution for the decay is small: 
r/o(i370)->pp-i-47r — 10 ± 10 MeV. (The error comes from varying Z between 1.6 and 1.88.) Ref. [I?} quotes 54 
MeV for the total 4-7T decay width. Since Ref. [48[ ascertains that about 26% of the total 4tt decay width originates 
from the pp decay channel, our result is consistent with these findings. 



V. CONCLUSIONS AND OUTLOOK 



We have presented a linear sigma model with vector mesons and global chiral invariance. The motivation for 
considering global invariance rather than the standard Sakurai model with local chiral invariance (with the exception 
of the vector-meson mass term) was that the latter fails to describe some important low-energy meson decay processes 
correctly, most notably the two-pion decay width of the p meson [l(|. This Lagrangian describes mesons as pure 
quarkonium states. As shown in Sec. IV. A, the resulting low-energy phenomenology is in general in good agreement 
with experimental data - with one exception: the model fails to correctly describe the /o(600) — ¥ irir decay width. 
This led us to conclude that /o(600) and cto(980) cannot be predominantly qq states. 

Assigning the scalar fields a and do of the model to the /o(1370) and ao(1450) resonances, respectively, improves the 
results for the decay widths considerably. We have obtained T^n^o^^ ~ 300-500 MeV for mj ( 1370 ) = 1200-1400 
MeV (see Fig. [HJ) - Thus, the scenario in which the scalar states above 1 GeV, /o(1370) and ao(1450), are considered to 
be (predominantly) qq states appears to be favoured over the assignment in which Jo (600) and do (980) are considered 
(predominantly) qq states. However, a more detailed study of this scenario is necessary, because a glueball state with 
the same quantum numbers mixes with the quarkonium states. This allows to include the experimentally well-known 
resonance /o(1500) into the study. 



19 



Of course, interpreting /o(1370) and ao(1450) as qq states leads to question about the nature of /o(600) and ao(980). 
Their presence is necessary for the correct description of tttt scattering lengths that differ from experiment for too 
large values of the isoscalar mass (see Sec. IV. A. 1). We distinguish two possibilities: (i) They can arise as (quasi- 
)molecular states. This is possible if the attraction in the irir and KK channels is large enough. In order to prove 
this, one should solve the corresponding Bcthc-Salpeter equation in the framework of Scenario II. In this case /o(600) 
and eto(980) can be classified as genuinely dynamically generated states and should not appear in the Lagrangian, see 
the discussion in Ref. [14j . If, however, the attraction is not sufficient to generate the two resonances /o(600) and 
ao(980) we are led to the alternative possibility that (ii) these two scalar states must be incorporated into the model 
as additional tetraquark states. In this case they shall appear from the very beginning in the Lagrangian and should 
not be considered as dynamically generated states. Of course, the isoscalar tetraquark, quarkonium, and glueball will 
mix to produce /o(600), /o(1370), and /o(1500), and the isovector tetraquark and quarkonium will mix to produce 
a (980) and <z (1450). 

The issue of restoration of chiral symmetry at nonzero temperature and density is one of the fundamental questions 
of modern hadron and nuclear physics, see, e.g., Refs. @,E3- Linear sigma models constitute an effective approach 
to study chiral symmetry restoration because they contain from the onset not only pseudoscalar and vector mesons, 
but also their chiral partners with which they become degenerate once the chiral symmetry has been restored. Once 
vacuum phenomenology is reasonably well reproduced within our model, we also plan to apply it to studies of chiral 
symmetry restoration at nonzero temperatures and densities. 

Another important check of the model is the description of the ALEPH data for the decay of the t lepton into two 
and three pions [49| . In this way, we will have a better constraint on the parameters of the model, e.g. the value for 
the a% mass. 

An extension of the model to Nj = 3 can be performed [5(|; with the exception of the strange quark condensate, 
no further free parameters will arise in this extension. However, much more data are available for the strange mesons, 
which constitute an important test for the validity of our approach. 

Acknowledgments 

The authors thank J. Schechter, F. Sannino, J.R. Pelaez, S. Gallas, and S. Striiber for valuable discussions. The 
work of D.P. was partially supported by the Foundation "Polytechnical Society". The work of D.P. and F.G. was 
partially supported by BMBF. The work of D.H.R. was supported by the ExtreMe Matter Institute EMMI. This work 
was (financially) supported by the Helmholtz International Center for FAIR within the framework of the LOEWE 
program launched by the State of Hesse. 

Appendix A: The full Lagrangian 

This is the final form of the Lagrangian (fT]) that is obtained after the shifts © and the renormalisation of the 
pseudoscalar wave functions; pP v = d^p u — d v p^\ a^" = d^a" — d^a^; (A)^ marks the third component of the vector 
A. Note that the term £4 contains the (axial-)vector four-point vertices [the terms ~ 54,5.6,7 in the Lagrangian ([T])]. 
We do not give the explicit form of £4 because it is not relevant for the results that are presented in this paper. 
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i (d^a + g 1 Zn-a t t+ giwZ 2 d^n ■ tt + giZrjf? + gi wZ 2 rj d^g) 2 



ml - c + 3 ( Ai + y 



a 2 



+ - {Zd^rr + gi Z^ x tt - 9l f?a - g lW Zd^r]a - 9l aa^ - g x wZad^^f 
+ i (Zd^g - gi af? - gxwZad^r} - g x 3? • a - giwZ ■ a f 



, , A2 
m - c + ( Ai + — 



rriQ + c + ( Ai + ^ 



Z 2 r? 2 



+ ^ [^a + ffip" x d + gi Zf?7r + 9l wZ 2 TTd^ + 9l Zg a? + 9l wZ 2 g d^f 



ml + c+ ( Ai + -A 2 



So - y [(™o + ^tt) 2 + Z 2 a 2 n 2 Z 2 (d ■ tt) 2 } 



1 I x , ^2 \ , 2 1 ->2 ^22 1 *72-;2\2 ( \ , ^2 \ , ,2 , -fl , ^2 2 , ^2=:2\ 



j Ai + yj (a 2 + a 2 + ZV + - ^Ai + yj 0a (a 2 + a 2 + Z 2 g 2 + Z 2 tt 2 ) 

- A 2 <M) • (™ + Z 2 gn) - 1 - dW) 2 + (c^ 1 ) 2 

- J - 0"^* + g 2 p f * x + g 2 a» x a1+ g 2 wZd^TT x d\ + g 2 wZa» x d L 'tt 
+ g 2 w 2 Z 2 (d f "7f ) x (d"ir)} 2 

+ - \ (**/r - ^/r) 2 + (f?r + (a?)- 

- i [0"a? - fl"^ + g 2 p" x d\ + g 2 wZp» x d v TT + g 2 dl{ x p" + g 2 wZ{d»Tr) x p"] 2 

- gj<j>d^ -[p^xZif- f?a - wZd d»g] 

- g\wZ <f> d^TT ■ [Zp>* x tt - /fob - wZd^g a ] + g\ (j>f lfl (a? ■ a + wZd^n ■ a ) 
+ g\wZ 4>d^g (d± ■ do + wZd^n ■ ao) 

+ g\ cfxj [(/f) 2 + IwZh^n + w 2 Z 2 {d»g) 2 ] 
+ g\(j}(T[{d>i) 2 + 2wZa lt , ■ d^TT + w 2 Z 2 (d»Tr) 2 ] 

- -^Z 2 {d,gd^ + d^-d^) 

z " l a 1 

+ eA^{(ao x c^a ) 3 + Z 2 {tt x d^) 3 - 4(p^ x p v ) 3 - A[{d lv + Zw8 v tt) x af] 3 

+ gi{2Z(f? + Zwd^g)(a x tt) 3 + Z{a + 0)[(a^ + ZwO^tt) x tt] 3 

+ g[a Q x (3? + Zwd^Tf)} 3 - al(a Q p^ + a 2 p^ 2 ) - ZV^V 1 + 7r 2 p^ 2 ) 

+ p^[( a J) 2 + (a 2 ) 2 + Z 2 (^) 2 + Z 2 ( 7 r 2 ) 2 ]} 

+ 4.g 2 {[pt + «L + ^ 2 ^ 2 (d^f + Zwa lu ■ 8 v tt]p^ + 2p v ■ (al + Zwd v Tt) x (af + ZwO^tt 3 ) 
-{p v -p^ + a lv ■ + Zwa lv ■ d^TT + Zwa^ ■ O u tt + Z 2 w 2 d u TT ■ ^ttV 3 

- [p^ ■ a lu + a^-p„ + Zwp^ ■ d v TT + Zwp v ■ d»w)af}} 

+ e lA^[(al) 2 + (a 2 ) 2 + Z 2 ^ 1 ) 2 + Z 2 (tt 2 ) 2 + 4(p^) 2 + 4(^ 2 ) 2 
+ 4(a{ v + Zwd^TT 1 ) 2 + 4(a 2 lv + Zwd v TT 2 ) 2 ] 

- 2e 2 A^A u [p^ p ul + p^ 2 p u2 + (af + Zwd^n 1 )^ + Zwd^ 1 ) 
+ (af + Zwd^TT 2 ){a v2 + Zwd v TT 2 )\ + £ hl23 + C 93 + U 
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+ 2Zw(f 1 ^i 1 + ai M • d»w)} 

(h\ h 2 h 3 \ 2 2 _z (hi h 2 h 3 \ 2 2 2 . 

+ (to - fc 3 )[(a + + ^V] • [f lft a1 + Zw(ai^ V + fi^Tf) + Z 2 w 2 (d^){d^)] 
+ (h 2 + h 3 )Z(do x 7r) • (w M a^ + Zwuj^d^ir) 
+ (h 2 - h 3 )Z{a x fr) • + Zwp^ri) 

+ h 3 Z[r)a - (a + ^)tt] • [p M x (a^ 1 + Zwd^Tv)] 

- y{(ao x p* 1 ) 2 - [Sb x (a? + Zu^tt)] 2 + Z 2 (tt x p* 1 ) 2 - Z 2 [vr x (a^ 1 + Zwd^n)} 2 } 



C g3 = -Ag 3 {d^[^ + f?fi +fr-(? + ^-a» + Zw(f?d v r) + f^r, + a? ■ d v n + a\ ■ 8^) 
+ Z 2 w 2 (d^ v d ,/ f] + d^n -d v n)} 

+ {dji„ + Zwd^rfil^fl + uj v ff( + ^ ■ a\ + j? ■ 3? + Zw{u»d v ri + u v d^ + ? ■ d v n + f? ■ 

+ d^p v • V + u V + /fa? + fffi + Zw(a1d v r] + ald^r) + f?d"n + ftd^n) + Z 2 w 2 (d»r)d v K + d v r}d fi 7r)] 

+ {d^ai v + Zwd^n) ■ [f?f? + ftp* + u^al + u v W{ + Zw{p^d v r] + fd^r) + u^n + uj v d^)]} 

+ Aeg 3 A^{u v [{F x jf) 3 + (3? x a^) 3 + Zw^n x a\) 3 + Zw(a» x d v T?) 3 + Z 2 w 2 (d^n x d v n) 3 ] 

+ {fiu + Zwd v r,)[{^ x d\) 3 + (3? x (?h + Zw{F x d v if) 3 + Zw{d^ x (?) 3 ]} 
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